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Abstract 

In this work, we derive the nuclear form factor for the spin-independent collision 
between the WIMPs and nucleus in terms of the relativistic mean field (RMF) 
theory. Comparison with the traditional form factors which are commonly used in 
literature is given and it is found that our results are slightly above that of the 2PF 
model by 4% to 8%, but deviate from the Helm form factor by 15% to 25% for 
the whole recoil energy spectrum of 0~ 100 keV. Moreover, taking Xe and Ge as 
examples, we show the dependence of the form factor on the recoil energy. 



1 Introduction 

Discovery of dark matter undoubtedly was one of the greatest scientific events of the 
20th century, then directly searching for dark matter and identifying it will be the most 
important task of this century. Commonly accepted point of view [Ij is that the main 
fraction of the dark matter in our universe which resides at the hallo of the galaxy is the 
cold dark matter, i.e. the weakly-interacting-massive-particles (WIMPs), and the most 
favorable candidate of the WIMPs is the neutralino [2], even though one cannot exclude 
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other possibilities at the present (Sj H]. What we expect is to "see" the dark matter 
particles in direct search experiments and find them at the LHC O |6[ [7] . The important 
step is to set experiments to directly search for such mysterious particles. Since the dark 
matter particles only participate in weak and gravitational interactions, the scattering 
cross section with normal matter is very small and makes the observation very difficult. 
Moreover, even though the dark matter particles are very heavy, its dispersive velocity is 
only 200 ~ 600 km/s [8], thus the available kinetic energy M(f /c)^/2 is of order of keV. 
This kinetic energy is not enough to make any inelastic nuclear reaction such as level 
transitions, therefore the only effects which are observable are the phenomena related to 
the nucleus recoil. This makes the measurement even more difficult in comparison with 
the neutrino experiments. 

The observation rates for spin-independent scattering can be written as [9] 



where Gp is the universal Fermi coupling constant, ao is the cross section at zero-recoil, 
rrir is the reduced mass of the WIMP and nucleus, finally F^(|q|) is the nuclear form 
factor. Generally speaking, the mass density distribution of the nucleus is proportional 
to the charge density or the nucleon number density, hence the form factor can also be 
accounted from the nucleon number density in the nucleus, i.e. the nuclear density. The 
form factor is the Fourier transformation of the nuclear density as 



where p{r) is the nuclear density. Here we assume that the nucleus is spherically symmetric 
and it is only a function of r. 

There are several ansatz for determining the form factor by assuming typical p(r) 
functions flU\ [TT| [T2| [13] . We will give a brief review of the more commonly adopted form 
factors in next section. 

The importance of the form factor is obvious. In fact, one of the the goals of particle- 
cosmology is to identify the dark matter. The strategy is to theoretically calculate cxo in 
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terms of various models for the WIMPs, which indeed determine how WIMPs interact 
with quarks in nucleon. The concrete models possess several model parameters [El [15] . 
Then one can compare the theoretical results with the data to extract information about 
the interaction and identity of the dark matter particle which is included in ctq. However, 
in Eq. (II]), ctq is convoluted with the nuclear form factor F{q) which is required to obey 

F(0) = 1. 

Unless one can more accurately determine F{q), extraction of useful information from 
data is impossible. Thus it is crucially important to determine the form factor. Since 
the recoil energy spectrum may span an energy range from a few hundreds eV to a few 
hundreds keV and the shape of the form factor may have different response to different 
recoil energies, so it is extremely important to have an accurate F{q) for the full energy 
range. Namely, an accurate determination of not only the values, but also the shape of 
the F{q) is necessary. 

In this work, we employ the relativistic mean field (RMF) theory to determine the 
nuclear form factors for Xe, Ge and several other nuclei which can be used as the materials 
for the dark matter detectors. Thus one can first choose any potential as the initial input, 
by iteration he eventually finds a stable form which makes all relevant equations to be 
self-consistent (see the text for details), the resultant form does not depend on the initial 
input which can be pre- assumed by our knowledge and intuition of physics. No doubt, if 
the initial potential is properly chosen, the iteration would converge quickly. Generally 
such a resultant form factor cannot be analytical, but only be presented numerically. 
With powerful computers, its application for extracting information of o"o from data is 
not affected at all. 

The article is organized as follows, after this introduction, we briefly review the forms 
of the nuclear form factors which are commonly adopted in literature, and then we derive 
necessary formulae for obtaining the form factor in terms of the RMF theory. In later 
section, we will present our results and make a comparison with the Helm and 2PF form 
factors in several figures. The last section is devoted to the discussion and analysis. 
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2 Nuclear form factor for the spin-independent cou- 
pling 



The nuclear form factor is the Fourier transformation of the nuclear density. In this 
paper, we numerically calculate the nuclear density by means of the RMF theory. Namely, 
we are going to derive the nucleon number density (nuclear density) in Eq. in terms 
of the RMF theory. 

2.1 Several types of nuclear form factors 

Here let us briefly review the form factors for the spin independent WIMP-nucleus 
interaction which were commonly adopted in literature. 

1. The simplest form factor is the exponential form [161 

F{E) = exp(-E/2Qo), (3) 

where E = /2Ma is the transferred energy from the WIMP to the nucleus and Qq = 
1.5/ (MRl). The authors of Ref . [T6l ITT] determined the nuclear radius Rq as: 

Ro ~ [0.9lA^/3 + 0.3]10-^3cm. 

2. The form of the two parameter Fermi (2PF) distribution is fl3[ ITS] : 

= eC'-^v" +r 

where po is equal to 2p(r) at r = c and the surface thickness t = 4a In 3. The parameters 
c and a for different nuclei have been determined by fitting the elastic electron scattering 
experiments and muonic atom spectroscopy [T9]. The Fourier transform of the 2PF density 
Eq. (jlj) cannot be given as an analytical expression, but can only be presented numerically. 
In literature [18] it is named as the 2PF form factor. 

3. The 'folding' charge distribution was given by Helm [IQJ as: 

p{r) = J po{r')pi{r -r')d^r', (5) 
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with certain normalization conditions. The nuclear charge density 



p[r) = po = SZe/AnR^, r <R 

is an approximately uniform distribution inside a cutoff radius R, and zero outside. The 
distribution pi(r) = (1/(2775^)^/^) exp(— r^/2s^) is introduced to take care of the soft edge 
of the nucleus where s is the surface thickness. Convoluting p(r) and pi(r), one deduces 
the Helm form factor [20J: 

F(q) = 3ji(gi?) ^_(g,)2/2 
qR 

^ sm{qR)-qRcos{qR) _(g,)2/2 ... 
" ^ (gi?)3 ^ ' 

with ji being the first-order spherical Bessel function. The Helm form factor is the form 
factor which is most commonly used in directly searching for dark matter [181 ED l22] . 
Lewin et al. [I3] proposed a set of parameters which are fixed by fitting the muonic atom 
spectroscopy [IH] as: 

R^ = c^ + IvrV-Ss^ (7) 
o 

c ~ 1.23A^/2 - 0.6 fm (8) 
s ~ 0.9 fm, a ~ 0.52 fm. (9) 

4. Sick [11] proposed to use a sum of the Gaussians (SOG) to express the nuclear 
charge density: 

where 7 is the Gaussians width. i?j and Qi are the position and charge fractions of the 
\—th Gaussian function with YliQi — 1- The three parameters are determined in the 
elastic electron scattering experiments [23]. The analytical form of the form factor is 
given as [18] : 

T-,, N _2 2/4^ Qi / D^ I 2_Rj^ sin(g_Rj) 
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5. In another model, the charge density is described as a sum of the Fourier-Bessel 
functions up to a certain cutoff radius R. The coefficients a,j of the Bessel expansions are 
directly extracted from the elastic electron scattering experiments. The charge density 
and form factor are given [121 [18] : 



2.2 Nuclear density in the RMF theory 

We employ the RMF theory to calculate the nuclear density in order to get the form 
factor for several nuclei which may be used as the materials of the dark matter detector. 
The RMF theory, which was first introduced by Walecka [23] , has been successfully applied 
in various fields of nuclear physics, especially the description of nuclear ground-state 
observables [23 ESI [23 [2B]- We adopt the RMF theory with nonlinear a and u terms, 
which are known to provide an excellent description of the ground states of finite nuclei 
including unstable nuclei [29| [30] . The Lagrangian for an ordinary nucleus in the RMF 
theory can be written as [29] : 

1 ^ ^ 1 

-^F^uF^" - ef^ATT^JeA^^^, (14) 
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with 



(15) 
(16) 
(17) 

where \I/Ar is the nucleon field with the mass M^r, a, cu^, are defined as meson fields 
with masses of m„,mi_j,mp. g„, g^^ and gp are respectively the couplings aNN, uNN, 





= dpUu 


- duUp,, 


Rflu 


= dpp^ 


- duPt,, 


F 


= dpA, 
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and pNN. g2, gs and C3 are coupling constants of the self-coupling terms for the a and 
CO mesons. is the electromagnetic (EM) field, r is the nucleon isospin operator and 
If. = {l + T-i)/2 is the projector and is the third component of the isospin Pauli matrices 
for nucleons. In the RMF theory, a, u, p and the EM field A are classical and treated as 
mean fields. In the Lagrangian, there are five independent fields: \1/Ar, a, u, p, and the EM 
field A. Employing the Euler-Lagrange equation, one can obtain five coupled equations 
corresponding to the five independent fields. Directly solving the equations is impossible 
because they are non-linear, thus one needs to introduce suitable approximations. These 
coupled equations can be solved in a standard way within the RMF approach [25]. 

We deal with static and spherically symmetric nuclei, and the meson fields depend 
only on radius, then the zero components u^, p^^ and A^ would make the dominant 
contributions. The equation of motion for the nucleon field is the Dirac equation where 
the classical fields a, u, p, and A stand as the effective potential for "^jy, 



{zj^d'' - [Mm + gMr)] - g.loUJ^r) - (7p7or3P°'(r) - e7o/cA°(r)}^^(r, t) = 0. (18) 



In the RMF theory the baryon field is still quantum. One can input a trial potential 
into Eq. (fTSl) and solve \l/iv, then substitute this solution into the equations for the 
meson and EM fields as the source terms to solve those fields. Iterating many times, the 
solutions would converge. The trial potential in principle can be freely chosen, but the 
convergence rate depends on the form of the trial potential. One can always choose a 
potential according to his knowledge on the nuclear structure or a reasonable intuition. 

We seek for solutions for \I^Ar(r, t) which have fixed energies. The nucleon field operator 
\1/Ar(r, t) can be written as a sum of positive and negative solutions ^^(r) and z/Q,(r) where 
a stands for the complete quantum numbers of an energy eigenstates and in the second 
quantization picture the wave function of a fermion can be written as 



Qa and h'^ are annihilation and creation operators of baryon and anti-baryon. The sub- 
script a runs over the full set of single-particle quantum numbers. For the normal nuclei, 
generally the negative solution z/a(r) can be ignored in the no sea approximation [25] . 





(19) 



a 
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The probability operator of the nucleon density at r is 

*jv(r,t)*jv(r,t). 



(20) 



The nuclear ground state can be expressed as 

>= (ntial)|o>, 

where |0 > is the vacuum, A is the number of nucleons in the nucleus and the product 
corresponds to creating A nucleons (protons and neutrons) which are the energy eigen- 
states. Sandwiching the probability operator Eq. fl20|) between the nuclear ground state, 
we obtain the nuclear density as 

ps(r) = ^4(r)«,(r). (21) 

a 

The input parameters are = 938 MeV, = 511.198 MeV, =783 MeV, 
= 770 MeV, = 10.0289, = 12.6139, gp = 4.6322, g2 = -7.2325, c/3 = 0.6183, 
and C3 = 71.3075, called as the TMl parameter set, which was determined by fitting 
some experimental data of finite nuclei in Ref . [29] . This parameter set is able to describe 
many ground state properties of finite nuclei all over the periodic table. We calculate the 
nuclear densities of several nuclei from light to heavy. In Fig. ([1]), we show the results of 
nuclear densities of ^^O, ^SNa, ^°Ca, ^^Ge, ^271^ i32xe, ^08pY^[_ 

Our next task is to perform the Fourier transformation of the RMF density numerically 
to gain the form factor. 

The kinematics for the collisions between WIMP and nucleus: the recoil energy is 

E = -mv^- TT-T^(l - COS^cm = " COS ^^m , 22 

2 {m + MAy Ma 
and the recoil three-momentum is 

|g1 = ^/2MAE, (23) 

^It is well known that the RMF theory is not a good working framework for even-odd nuclei, such as 
^■^Na and ^^^j^ since those elements are possibly chosen as the dark matter detector materials, we use 
the RMF approach to evaluate their form factors as a reference for further studies. We will investigate 
their form factors in a more proper framework, if such elements are indeed chosen to be the detector 
materials. 
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Figure 1: Nuclear densities of ^^O, ^^Na, ^°Ca, "^Ge, ^^"l, ^^^Xe, ^°^Pb, obtained in terms 
of the RMF theory. 



where m and Ma are the masses of the WIMP and nucleus respectively, is the reduced 
mass, V is the velocity of the WIMP and 6cm is the scattering angle in the center of mass 
frame |3lj- In the lab frame, the recoil momentum is exactly the momentum transferred 
from the incident WIMP to the nucleus which is supposed to be at rest before impact of 
the WIMP. 

Taking Ge as an example, in Fig. ([2]) we show the dependence of the RMF form factor 
on the recoil momentum. For a comparison we also draw the curves corresponding to the 
Helm and 2PF form factors. The relevant parameters of the Helm and 2PF are taken 
from Refs. [iHl |19]. As can be seen from the Fig. ([2]), the results of RMF are 2% to 8% 
above that of 2PF form factors at most of the energy region from 2 keV (momentum q ~ 
0.1 fm~^) to 200 keV (q ~ 0.9 fm~^), and lower by less than 1% at the range of g < 0.1 
fm~^. But at all the recoil energy regions, the RMF results are above the Helm form 
factor by about 15% to 25%. At 1 fm~^, all of them attenuate to zero, or very close to 
zero. 
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Figure 2: RMF, 2PF and Helm form factors for ^^Ge for the recoil energy range from 10 
keV (q ~ 0.19 fm'^ ) to 100 keV (q ~ 0.61 fm^^ ). 




Figure 3: The RMF form factors for various nuclei corresponding to the nuclear densities 
shown in figure 1. The left and right pictures are the dependence of the factors as functions 
on respectively the recoil energy and the recoil momentum. 
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We plot the dependence of the RMF form factors on the recoil energy and recoil 
momentum (for a clearer illustration) for various nuclei in Fig. ([3]). It is shown that the 
form factor generally drops faster as the nucleus mass increases at higher recoil energies. 
Among the commonly adopted detector materials, Fxe(E)(dash dot line) and F/(E) (short 
dash line) have similar behaviors. They drop off more rapidly at the recoil energy range 
from 10 keV to 100 keV and Fxe(E) tends to zero at 60-70 keV. It significantly influences 
the detection sensitivity of the Xenon detector. At the recoil energy larger than 100 keV, 
F7-(E) and FGe(E)(dot line) are too small to perform effective detections. In the case of 
Na and Ca, the form factors decrease more slowly in the whole energy spectrum. 



3 Discussion and Conclusions 



Though the RMF form factor is not analytical, its advantage is obvious. We are based 
on the theory which has been proved to be successful for dealing with nuclear problems, 
especially determining the nuclear densities for spherically or near-spherically symmetric 
nuclei. In the RMF theory, by resolving the Dirac equation to get nucleon wave function 
and with it the nuclear density can be obtained. Since this method does not depend on 
the initial form of the potential, even though we need to choose one for the first input, it 
has obvious advantage. The iteration would smear the difference between different inputs 
and make the final solution to converge to a should-be reliable one. Therefore we believe 
that the RMF form factor would be a more accurate form factor for the calculation of 
dark matter detection, and moreover, it also determines the energy-sensitivity curve of 
the material adopted for the detector. 

The results can be utilized as a reference for designing dark matter detectors and make 
accurate extraction of information about the identity of the WIMP and its interaction 
with nucleon, or even with the ingredients inside nucleons (say, quarks) from data j32l [33t 
[32 Eg [36]. 

For readers' convenience, besides the figures we presented above, we make a table 
where some numbers are explicitly listed. As can be seen in Table 1, at the recoil energy 
ranges, the form factors of the three models have ceratin differences, the RMF and 2PF 
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form factors differ from each other by only 4%-8%, but differ by at least 15% to 25% from 
the Helm form factor. The Fig. ([3]) tells us from 10 keV to 70 keV, Fxe{E) obviously 
declines from FceiE) and Fj{E). 



Table 1: the three form factors at 10 keV, 50 keV, 100 keV , for different nucleus. 



E 


\nE)? 


016 


Na23 


Ca40 


Ge72 


1127 


Xel32 


Pb208 


10 keV 


RMF 


0.982 


0.968 


0.926 


0.818 


0.596 


0.576 


0.293 


2pF 


0.981 


0.965 


0.920 


0.802 


0.605 


0.579 


0.305 


Helm 


0.908 


0.882 


0.825 


0.710 


0.507 


0.489 


0.251 


50 keV 


RMF 


0.913 


0.852 


0.678 


0.351 


0.047 


0.035 





2pF 


0.907 


0.838 


0.658 


0.324 


0.052 


0.042 





Helm 


0.777 


0.706 


0.545 


0.278 


0.040 


0.031 





100 keV 


RMF 


0.834 


0.704 


0.454 


0.105 











2pF 


0.809 


0.700 


0.429 


0.098 











Helm 


0.671 


0.567 


0.346 


0.083 












Here, We assume the target nuclei discussed in this work to be spherical or near- 
spherical, however realistic nuclei may be slightly deformed. Such deformations possibly 
affect the shape of the form factors and we are going to discuss the effects in our next work. 
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